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Abstract 

A set of analytical benchmarks for tracking programs is required for precision storage ring experiments. To determine 
the accuracy of precision tracking programs in electric and magnetic rings, a variety of analytical estimates of particle 
and spin dynamics in the rings were developed and compared to the numerical results of tracking simulations. Initial 
discrepancies in the comparisons indicated the need for improvement of several of the analytical estimates. As an 
example, we found that the fourth-order Runge-Kutta/Predictor-Corrector method was slow but accurate, and that it 
passed all the benchmarks it was tested against, often to the sub-part per billion level. Thus, high precision analytical 
estimates and tracking programs based on fourth-order Runge-Kutta/Predictor-Corrector integration can be used to 
benchmark faster tracking programs for accuracy. 

Keywords: Analytical Benchmarking, Precision Particle Tracking, Electric and Magnetic Storage Rings, Runge-Kutta, 
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1. Introduction 

Analytical estimates for particle dynamics in electric 
and magnetic rings with and without focusing have been 
given in a variety of papers and notes. We have used these 
high precision estimates to provide benchmarks to test the 
accuracy of precision particle tracking programs. The term 
“focusing” in this paper is used to denote “weak vertical fo¬ 
cusing” unless otherwise indicated. Horizontal focusing is 
defined by the vertical focusing plus the geometry of the 
ring, always conforming with Maxwell’s equations. 

These benchmarks include: 

• Pitch correction (TJ [2] to particle precession frequency 
in a uniform B-field with and without focusing. 

• Vertical oscillations and energy oscillations in a uni¬ 
form B-field with no focusing, electric focusing, and 
magnetic focusing. 

• Radial and vertical oscillations and energy oscillations 
in an all-electric ring with and without weak focusing. 

• Synchrotron oscillations and momentum storage with 
a radio frequency cavity (RF) in a uniform B-field. 

• An EDM signal and systematic error with an RF Wien 
Filter in a magnetic ring. 

In the analytical estimates that follow, we define 70 
as the Lorentz factor of the design particle in the ring. 


The vertical pitch angle 6 y of a particle is defined such 
that Oy = Pz/Pq where (/3 p , /?#, /3 Z ) = v/c in cylindri¬ 
cal coordinates. The field focusing index is n, with n = 
— (dB/Bo)/(dr/ro) a number with range 0 < n < 1. 

These analytical estimates provide a means of bench¬ 
marking particle tracking programs. A precision track¬ 
ing method that successfully passes all benchmarks can 
provide a baseline to benchmark faster programs. We 
study the developed benchmarks in the context of Runge- 
Kutta/Predictor-Corrector integration. The validity of the 
developed benchmarks extends to any precision tracking 
program. 


2. Motivation 

Precision experiments such as the Muon (g— 2 ) and Stor¬ 
age Ring Electric Dipole Moment (EDM) experiments|3, 
HE] require measurements of sub-part per million (ppm) 
accuracy. In the case of a proton or deuteron Storage 
Ring EDM experiment, a tracking program of extraordi¬ 
nary precision is required to estimate the spin coherence 
time of the particle distribution and various lattice param¬ 
eters, as well as to estimate the values of systematic errors 
associated with the experiment. Many commonly used 
beam and spin dynamics programs ignore, or erroneously 
account for, second and higher-order effects. Tracking in 
an electric storage ring poses the additional challenge of 
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conforming with total-energy conservation while account¬ 
ing for higher-order effects. 

Numerical integration with a sufficiently small step size 
allowed to run for a sufficiently long time may reproduce 
the analytical result with high accuracy. A tracking pro¬ 
gram to be used for estimates and investigations in preci¬ 
sion experiments must be optimized to be as fast and ac¬ 
curate as possible. This calls for a well-tested and robust 
procedure to benchmark the accuracy of tracking programs 
in situations relevant to the experiments. 

We summarize and derive analytical solutions to the 
equations of motion of a particle in various electric and 
magnetic rings. In several cases, comparison of the an¬ 
alytical estimate with precision tracking results identi¬ 
fied discrepancies and indicated the need to improve the 
estimates. In this way it was determined that the to¬ 
tal correction due to vertical particle oscillations, the so- 
called pitch effect, can be significantly reduced[6]. These 
analytical estimates provide individual benchmarks for 
tracking programs. A program well-benchmarked against 
these estimates would provide additional, more flexible 
means to benchmark faster programs in different con¬ 
ditions. We benchmarked a program based on Runge- 
Kutta/Predictor-Corrector method [7] against the devel¬ 
oped analytical estimates. 

Runge-Kutta/Predictor-Corrector integration should 
reproduce the analytical estimates to sub-ppm accuracy 
on a time scale on the order of hours, in order to be 
a feasible candidate for benchmarking faster programs. 
Together, the analytical estimates and a program based 
on Runge-Kutta/Predictor-Corrector integration provide 
a powerful tool for benchmarking precision tracking pro¬ 
grams for Muon (g — 2), Storage Ring EDM, and other pre¬ 
cision physics experiments requiring high precision beam 
and spin dynamics simulation. 


3. Precision Tracking 


For a particle of mass m and charge e, there are two dif¬ 
ferential equations that govern particle and spin dynamics. 
For particle velocity (3 and rest spin s in external fields, 
the equations are[8 : 


d£ 

dt 


m^fc . 


E + c(3 xB-j30-E) 


and the T-BMT equation, with an anomalous magnei 
moment a of the particle: 
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The Predictor-Corrector integration was used with a 
step size of 1 — lOps to numerically solve the two differ¬ 
ential equations with the corresponding initial conditions. 


Although the method is relatively slow, it is very simple 
and accurate. This method uses the Runge-Kutta method 
to start the integration process and we refer to it as Runge- 
Kutta/Predictor-Corrector in this document. 


4. Magnetic Ring 

A magnetic ring consists of a uniform magnetic field R, 
taken to be in the vertical direction. The correction C to 
the precession frequency due to a vertical pitch is defined 
by cJm = cj a (1 — C), where uo a is the (g — 2) frequency [9] 
for a particle with anomalous magnetic moment a. The 
predicted correction is|2|: 

C = \°0 U - ( w a + 2a 7 2 Wp)/7 2 (^a - y)}. (3) 

with uj p = 27r/ p , where f p is the vertical (pitch) oscillation 
frequency. 


4-1- No Focusing 

When there is no focusing or when uo p cj a , the cor¬ 
rection from Equation [3] becomes: 

C = \p 2 el (4) 

where for linear oscillations, ( 0 y) = where Oq is the 
maximum pitch angle of the particle trajectory. 

For a particle with (3 = 0.972 and a constant l.Omrad 
vertical pitch as shown in Figure [l] the simulated correc¬ 
tion to the (g — 2) precession frequency of 0.2361ppm is in 
very good agreement with the analytically predicted value 
of 0.2363ppm using Equation |4j 

Checking over several values of 0 y confirms that the an¬ 
alytic expression and the pitch correction in the tracking 
simulation agree for small Q y , as expected. 


4-2. Weak Magnetic Focusing 

When there is magnetic focusing and when oo p u; a , 
the correction from Equation [3] becomes: 

c = 4 O 0 O- + 2 a). (5) 

The analytical estimate [TO] for the average particle ra¬ 
dial deviation from the ideal orbit with radius ro, with 
weak magnetic focusing index n, takes the form: 




1 


1 — n 




(6) 


for a vertical pitch frequency significantly greater than the 
(g — 2) precession frequency of the particle, where a p is the 
momentum compaction factor. 

Equation |6l predicts an average radial deviation (Ar/ro) 
of —5 x 10 ' using 6 o = lmrad and a field index n = 
0.01, consistent with the tracking results shown in Figure 
[5] to sub-part per billion (ppb) level. The dependence of 
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(Ar/ro) on the field index is shown to hold over a range 
of n values in Figure [3j 

In a continuous storage ring with weak focusing, field 
strength Bp, and ring radius rp, the vertical and horizontal 
magnetic field components around the ideal trajectory can 
be expressed to second-order in the vertical position y as: 

B x (x,y) =-n—y (7) 

r o 


T 3 ( ^ d B 0 Bq y 2 

B y {x, y) = B 0 — n—x + n — —, (8) 

r 0 r 0 2 r 0 

where the nonlinearity arises from the application of 
Maxwell’s equations in cylindrical coordinates. The hor¬ 
izontal and vertical tunes are given by v x = \/l — n and 
v y — yjn respectively. 

We make use of the relations 10 from Equation [6j 


80 = 1 #0 

rn P 2 1-n 2 


(9) 


and #o = Vo/ (ro/y/n), where yo is the maximum vertical 
excursion. From this we see that: 
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By considering the time-averaged relative B-field 
change, calling Cb the modification due to the differ¬ 
ent B-field encountered by the particle, we find that the 
correction C* = C — Cb to the (g — 2) precession fre¬ 
quency in a magnetic storage ring with weak focusing, i.e. 
l dm = u a (l — C*), is given by the expression: 


C* 



2 ’ 


( 11 ) 


where we see that several terms of the inhomogeneous 
B-field correction and the correction in Equation [5] can¬ 
cel, leaving a small correction. The necessity of including 
the second-order inhomogeneous magnetic field contribu¬ 
tions was overlooked by previous authors. Investigations 
with precision tracking identified the discrepancy and mo¬ 
tivated the improvement of the analytical estimate. 

This holds for a vertical pitch frequency much greater 
than the (g — 2) precession frequency of the particle, which 
for a weak focusing ring means y/n ay. Equation 11 


implies that the pitch effect can, in principle, be made to 
vanish for n = a/(l + a), but the condition y/n ay 
makes it rather difficult to achieve. To test the tracking 
program we introduced a particle with 10 times the muon 
mass, with same a value as the muon, stored in a ring 
radius of 7.112m. The program indeed showed that the 
pitch correction vanishes with an uncertainty at the part 
per billion (ppb) level when n = a/(l + a) was used. 

For realistic muon parameters, the observed (g — 2) fre¬ 
quency is off from its correct value by +0.109ppm, for a 
vertical maximum pitch angle Oq = lmrad, and n = 0.18, 


consistent with the offset shown in Figure [4] to sub-ppb 
level. 

A resonance of the pitch effect correction occurs when 
the vertical pitch frequency uj p is equal to the (g — 2) 
precession of the stored particles uj a . The correction C 
approaches Equation [4] for uj p u: a and Equation [5] for 
uj p uj a . The full range of pitch corrections C over a 
range of index values is shown in Figure [5j When all the 
fields are taken properly into account, as shown above, the 
tracking results reproduce the same curve to sub-ppb level 
for A Ld a /uj a . 

A comparison of the frequency shift predicted by Equa¬ 
tion [TT] and the results from tracking is given in Table [T] 
The analytical estimates of the pitch correction and the 
tracking results are in very good agreement, better than 
ppb level. This level of precision is adequate for the Muon 
g — 2 experiments currently underway El HI], both aiming 
for better than O.lppm total systematic error. 

4-3. Weak Electric Focusing 

In the case of electric focusing in a uniform magnetic 
ring, the expected precession frequency correction due to 
the pitch effect is [2]: 

c = \e\]{ p 2 - («wy )/u - y)}. (12) 

and 

^ = + a )i (13) 

for uo p uj a and for the particle at the magic 

momentum [4j [5j such that an electric field does not affect 
the (g — 2) precession. 

For a maximum vertical pitch of So = 0.5mrad, the ana¬ 
lytically estimated pitch correction of 0.0624ppm for magic 
momentum muons is very close to the result from tracking, 
shown in Figure [6j consistent to sub-ppb level. 

The results from tracking match the predicted value to 
sub-ppb accuracy. Thus, we conclude that the analytical 
estimates and the Runge-Kutta/Predictor-Corrector inte¬ 
gration method have passed the magnetic ring tracking 
benchmarks. 

4-4- Radio Frequency Cavity 

The synchrotron oscillation frequency f s of a particle in 
a uniform B-field with an RF is: 


fs = Qsfc, (14) 

where f c is the cyclotron frequency and Q s is the syn¬ 
chrotron tune, which satisfies: 

2 = eVprich 
s 2lTCpoP 2 ' 

In the expression above, e is the elementary charge, Vo 
is the voltage of the RF cavity, h is the harmonic of the 
RF cavity used, and po is the ideal momentum. The value 
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of Tj c , the so-called slip factor, is determined from the ex¬ 
pression: 


Tjc — CXp 
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l l 

l-n 7 g' 


(16) 


Using a particle with charge e, po = 3.094GeV/c, and 
7 o = 29.3, and using a 20cm RF cavity, with Vo = 100kV 
and harmonic h = 1 , the predicted synchrotron frequency 
with n = 0.18 is f s = 16.8kHz and with no vertical fo¬ 
cusing it is f s = 15.2kHz. Comparing these calculations 
with the results of tracking in Figure [7] shows close agree¬ 
ment between the tracking simulation and the estimation, 
verified at the 0 . 1 % level. 

The maximum momentum storage range of the RF 
cavity 1 12 ] is given by the expression: 


A p 
Po 


\ _ 1 

1 2 eV 0 

) V 
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TThrjcPcpo 


(17) 


around the ideal particle with momentum po. 

Using the above RF parameters and particle values, the 
maximum stored momentum is estimated from Equation 


The RF phase diagram shown in Figure [ 8 ] is consistent 
with the estimated value and illustrates the momentum 
storage region of the configuration. 

The analytical estimates for the maximum stored mo¬ 
mentum and the synchrotron frequency were matched 
closely by the Runge-Kutta/Predictor-Corrector tracking 
method. Thus, we see agreement to the desired accuracy 
between the analytical estimates and results of Runge- 
Kutta/Predictor-Corrector integration in the case of a 
magnetic ring. 


17 to be (Ap/po) max = 0.00454 with no vertical focusing. 


5. Electrostatic Ring 


In cylindrical coordinates, the electric field with an index 
m power law dependence on radius at y = 0 is: 

r l+m 

£M) = £o^f, (18) 


where y is the vertical direction and r is in the radial di¬ 
rection. 

In a uniform all-electric ring, we have found that the 
y —y and rotational symmetries allow the radial and 
vertical electric field components to be found exactly: 


E y (r 1 y) = E 0 r ¥ mmy 


pl-\-Tn 


_ , m „ m3 7 % , . 
\F\ (l + —, 1 + —; —; — ——) (19) 
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for all m > 0 , where 2^1 is the ordinary hyper geomet¬ 
ric function. The field index is n = m + 1, and m = 0 
corresponds to cylindrical plates with no vertical focusing, 
m = 1 corresponds to spherical plates, and so on. The fo¬ 
cusing m value depends on the choice of electrode profile. 


For m = 0, the electric field is that of a uniform cylin¬ 
drical capacitor. The fields were taken to fifth-order in y/r 
when implemented in the tracking program. The expan¬ 
sion to second-order is shown below: 


E r (r,y) = Eq 7 ^- 


i -V n2 - i ^ + ®w 


E y (r,y) = E 0 ^ 


(n — !)- + €?( ^ 


( 21 ) 

( 22 ) 


The contributions of the higher-order electric field terms 
were found to be negligible for tracking. The second-order 
term is significant for the analytical estimates. The fields 
given in Equations [19] and [20] describe the field configura¬ 
tion considered for an electric ring. 

In an all-electric ring, the kinetic energy changes with 
the radial position, which provides additional horizontal 
focusing. The horizontal and vertical tunes [13 14j are 
given by v x = \/l — ra + l/ 7 2 and v y = ^Jm respectively, 
in an electric ring with weak focusing. 


5.1. No Focusing, including an RF cavity 

With no focusing in the ring, we have n = 1 and thus 
m = 0, corresponding to concentric cylindrical plates. We 
also include an RF-cavity, which fixes the particle revolu¬ 
tion frequency. Y. Orlov [T5l [16] and I. Koop 17 solved the 
orbital motion for an electrostatic field with no focusing. 
In this case, the estimates for the average values of A 7/70 
and Ar/ro take the following form: 
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(24) 


Note that these values depend only on the particle ideal 
Lorentz factor 70 and the pitch angle, not on the ring 
radius, plate spacing, or electric field strength. 

The precision tracking results for the two parameters 


and the predicted values of Equation [23] and Equation [24 
are shown in Figure [9] We see close agreement between the 
expected value and the values calculated through tracking. 
Incidentally we found from tracking that), without includ¬ 


ing an RF-cavity, 


/ A r 

\ ro 


from Equation 


24 


becomes: 


Ar\ , 2 . 7o 

r 0 / {y) il + r 


(25) 


whereas Equation [23] remains the same. 


5.2. Weak Electric Focusing, including an RF cavity 

With weak focusing such that 0<m<l, the param¬ 
eters analytically estimated by Y. Orlov [15j [16] are given 
by Equations [26] and [27] below: 
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which hold for times much larger than the period of ver¬ 
tical oscillations. Note that these values depend only on 
the pitch angle and not on the ring geometry, ideal 70 or 
the field focusing index. 

There is an apparent gap between Equations [23] and 


transition between focusing and no focusing can exist since 
the latter formulas hold only for averages over times much 
larger than the period of vertical oscillations [16]. 

Figure [10] and Figure [Tl] shows the comparison of the 
precision tracking results with the analytical estimates for 
(A 7 / 70 ) and (Ar/ro), respectively. (There is a vertical 
spread of less than one part per billion for (A 7 / 70 ) and 
less than 0.1 ppm for (Ar/ro), which we assign as the error 
of the method.) 

The results from tracking in Figures [9]jTT] match the pre¬ 
dicted (A 7 / 70 ) and (Ar/ro) values to sub-ppm accuracy. 
From these comparisons, we see again that the analyt¬ 
ical estimates and the Runge-Kutta/Predictor-Corrector 
method pass the electric field benchmarks. 


24 and Equations and pH in the limit as m -A 0. The 


6. Radio Frequency Wien Filter 


A radio frequency Wien Filter (WF) is a velocity- 
dependent charged particle filtering device. A WF can 
be used in a storage ring to measure a particle’s electric 
dipole moment (EDM). The analytical estimatesp~ 8 ] of the 
EDM signal and systematic error for a particle of charge 
e, mass m, and anomalous magnetic moment a in electric 
and magnetic fields are: 


/ dsy\ eby (1 + a) e(— Er + /3By) 

V dt Jedm VAmC T 2 Si ° mCUJafl 


(28) 
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where slo is the peak longitudinal spin magnitude, Er 
is the radial electric field strength, and By is the verti¬ 
cal magnetic field strength. The EDM is proportional to 
7 , with r] playing the same role for the EDM as the g- 
factor plays for the magnetic dipole moment. The radio 
frequency WF ideally produces a vertical magnetic field at 
the g — 2 frequency by = byo cos uo a ot and a radial electric 
field e R = e R o cos u: a ot with the condition[18] e R o = /3byo. 
However, if the WF is misaligned by an angle 9 with re¬ 
spect to the vertical, then a radial B-field will also be 
present b R = b R o cosco a o£, with b R 0 = by 0 sin 6 inducing 
a systematic error given by Equation [29] 

A comparison between the analytical estimates and the 
tracking results for the deuteron case and the proton case 
are given in Tableland Table [3J respectively. 


The Wien Filter provides another benchmark for testing 
the accuracy of the analytical estimates and the Runge- 
Kutta/Predictor-Corrector tracking method. We again see 
very good agreement between the analytically predicted 
values and those calculated by tracking. 


7. Conclusion 

We have determined an array of analytical estimates for 
benchmarking tracking programs for precision storage ring 
experiments. The benchmarks form a robust test for elec¬ 
tric and magnetic rings, with and without focusing, as well 
as RF cavities and Wien filters. The inclusion of a high- 
order contribution to the weak magnetic focusing estimate 
resulted in the discovery of a method to reduce or eliminate 
the pitch effect, which was overlooked by other authors. 
Together these analytical estimates give a powerful tool 
to benchmark programs for studying particle motion and 
spin dynamics in a variety of storage ring configurations. 

The Runge-Kutta/Predictor-Corrector integration ap¬ 
pears to be an accurate, albeit slow, tool for precision 
tracking. It has passed all the benchmarks that it was 
tested against, often to the part per billion level. The 
tracking program was able to successfully simulate par¬ 
ticle dynamics in electric and magnetic rings with and 
without weak focusing, in agreement with analytical es¬ 
timation. Consequently, we conclude that the Runge- 
Kutta/Predictor-Corrector method can be used to bench¬ 
mark faster tracking programs. Other tracking programs 
are also capable of providing high accuracy at least under 
certain conditions, e.g., pa mi. What we have presented 
here is an array of high precision analytical estimates that 
can be used to benchmark candidate simulation programs. 
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Figure 1: The particle path in Cartesian coordinates in a uniform 
B-field with pitch angle 6 y = l.Omrad, for a ring with a 5m radius. 


Radial deviation vs. Time 



Figure 2: The particle deviation from the ideal radial position over 
time, modulo 50/iS. The simulations used a maximum pitch angle of 
6q = lmrad and magnetic focusing with field index n = 0.01. 



Figure 3: The average (Ar/ro) versus the field focusing index n. 
The solid line represents the predicted values while the points are 
the results of tracking. The simulation used a maximum pitch angle 
of 6> 0 = l.Omrad. 
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Figure 4: The pitch correction to the (g — 2) frequency with angle 
6 q = l.Omrad, 7 = 29.3, and n = 0.18, is 0.109ppm, consistent with 
Equation |11| to sub-ppb level. 


Figure 6 : The (g — 2) relative difference estimated from tracking using 
electric focusing in a magnetic ring with a maximum pitch angle of 
60 = 0.5mrad. To obtain the correct g — 2 frequency, the correction 
0.25 X 0q needs to be added to the observed frequency. Again, the 
tracking results are consistent to sub-ppb level with the predictions. 



Figure 5: Parameter C from Equation [3] for the pitch correction to 
the (g — 2) frequency over a range of n values using 6 y = 1.5mrad and 
/3 = 0.94. A resonance occurs when u p = uj a , at n « 1.25 X 10 —5 , as 
expected. 


Table 1: Comparison of the frequency shift estimated using Equation 
|ll| and the tracking results. The tracking results assume a muon with 
7 = 29.3, stored in a magnetic ring with magnetic focusing and a 
radius 7*0 = 7.112m. The vertical angle used is 0 y = 0.5mrad. The 
observed g — 2 frequency is shifted higher by the small factors given 
below depending on the n -value used. 


n 

Estimation (ppb) 

Tracking (ppb) 

0.01 

1.1 

1.0 

0.02 

2.4 

2.4 

0.03 

3.7 

3.6 

0.05 

6.4 

6.4 

0.08 

10.7 

10.8 

0.10 

13.7 

13.7 

0.137 

19.7 

19.9 

0.237 

38.7 

38.8 




Figure 7: The synchrotron oscillations of a particle in a uniform 
B-field with an RF with (top) n = 0.18 and (bottom) no vertical 
focusing. The solid line represents the estimated oscillations; the 
points are the results of tracking. 
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Figure 8 : A phase diagram of particle momenta in a uniform B-field 
with an RF and no vertical focusing. The closed energy oscillations 
around the synchronous particle define the storage region of the RF. 
The unstored particle is close to the boundary of the storage region. 
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Figure 9: The average (A 7 / 70 ) (on top) and average (Ar/ro) (on 
bottom) versus the ideal 70 of the proton. The solid lines represent 
the predicted values; the points are the results of tracking. The 
simulation used a pitch angle of 6 y = l.Omrad, an RF-cavity and no 
vertical focusing. 
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Figure 10 : The average (A 7 / 70 ) versus the ideal 70 of the proton over 
a variety of focusing n values. The solid line represents the predicted 
values; the points are the results of tracking. The simulation used 
0q = l.Omrad, an RF-cavity and vertical focusing. 
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Figure 11 : The average (Ar/ro) versus the ideal 70 of the proton over 
a variety of focusing n values. The solid line represents the predicted 
values; the points are the results of tracking. The simulation used 
#0 = l.Omrad, an RF-cavity and vertical focusing. 


Table 2: Radio Frequency Wien Filter: Comparison between analyti¬ 
cal estimates and tracking results for the deuteron case, in rad/s. Mo¬ 
mentum p is in GeV/c. The EDM value is assumed to be 10 _ 18 e-cm, 
while for the (systematic) error, the misalignment angle is assumed 
to be O.lmrad. 


p 

EDM 

Tracking 

EDM 

Analytic 

Error 

Analytic 

Error 

Tracking 

0.7 

-1.00 

-1.00 

0.41 

0.41 

1.4 

-0.74 

-0.73 

0.175 

0.17 

2.1 

-0.50 

-0.51 

0.096 

0.097 

2.8 

-0.36 

-0.35 

0.063 

0.06 
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Table 3: Radio Frequency Wien Filter: Comparison between analyt¬ 
ical estimates and tracking results for the proton case, in rad/s. Mo¬ 
mentum p is in GeV/c. The EDM value is assumed to be 10 —18 e-cm, 
while for the (systematic) error the misalignment angle is assumed 
to be O.lmrad. 


p 

EDM 

Tracking 

EDM 

Analytic 

Error 

Analytic 

Error 

Tracking 

0.7 

0.357 

0.357 

1.135 

1.137 

1.4 

0.174 

0.172 

0.393 

0.396 

2.1 

0.0934 

0.093 

0.192 

0.195 

2.8 

0.0566 

0.0563 

0.1135 

0.1135 


9 



